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We calculate the one- and two-loop corrections of order a (Z a) 6 and a 2 (Z a) 6 respectively, to the Lamb 
shift in hydrogen-like systems using the formalism of nonrelativistic quantum electrodynamics. We obtain gen- 
eral results valid for all hydrogenic states with nonvanishing orbital angular momentum and for the normalized 
difference of S-states. These results involve the expectation value of local effective operators and relativistic 
corrections to Bethe logarithms. The one-loop correction is in agreement with previous calculations for the 
particular cases of S, P, and D states. The two-loop correction in the order a 2 [Z a) 6 includes the pure two- 
loop self-energy and all diagrams with closed fermion loops. The obtained results allow one to obtain improved 
theoretical predictions for all excited hydrogenic states. 
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I. INTRODUCTION 

The precise calculation of the electron self-energy contri- 
bution to energy levels of hydrogen-like systems is a long- 
standing problem in bound-state quantum electrodynamics. 
The widely used direct numerical approach 01 Si is based 
on a partial-wave decomposition of the Dirac-Coulomb prop- 
agator, which corresponds to the exact all-order treatment of 
the electron-nucleus interaction. The one-loop corrections 
have already been calculated to a high numerical precision for 
a wide range of nuclear charge numbers Z (including the case 
of atomic hydrogen Z — 1), whereas the two-loop correction 
has been obtained only for Z > 10 with limited numerical 
precision. The analytic method is based on an expansion in 
powers of Z a and a subsequent analytic or semianalytic in- 
tegration. The two approaches are complementary. In prac- 
tice, the numerical method has primarily been used for sys- 
tems with a high nuclear charge number, whereas the analytic 
method usually provides more accurate predictions for low-Z 
systems. 

Here, we present a unified analytic derivation of the one- 
and two-loop binding corrections of order (a/V) (Z a) 6 m c 2 
and (a/n) 2 (Z a) 6 mc 2 respectively, for arbitrary bound 
states of hydrogen-like system using the formalism of dimen- 
sionally regularized nonrelativistic quantum electrodynamics 
(NRQED). This method allows for a natural separation of dif- 
ferent energy scales, (i) the electron mass and (ii) the binding 
energy, using only one regularization parameter: the dimen- 
sion d of the coordinate space. This leads to a straightforward 
derivation of radiative corrections in terms of expectation val- 
ues of some effective operators and the Bethe logarithms. The 
calculation of these operators is the main task of this work, 
and we obtain them from standard electromagnetic form fac- 
tors and the low-energy limit of the two-photon exchange scat- 
tering amplitude. 

This paper is organized as follows: In Sec.||lJ dimension- 



ally regularized NRQED is outlined. In Sec. Ill, the one- 
loop self-energy is derived by splitting the calculation into 
low- (Sec. |nn|}, middle- (Sec. lIIICl . and high-energy parts. 
The general one-loop result is presented in Sec. 1111 Dl and the 
evaluation for D, P and S states in Secs. liTOIillTl and lilTGl 
respectively. The two-loop correction is separated into four 
different gauge-invariant sets of diagrams, see Figs.^— E]be- 
low. These are subsequently investigated in Sees. IIV1 — IVIll 
Results are summarized in Sec. lIXI Moreover, in Appendix C 
we present the calculation of an additional two-loop logarith- 
mic contribution to the ground state which was omitted in the 
previous work |4j] . 



II. DIMENSIONALITY REGULARIZED NRQED 

As is customary in dimensionally regularized QED, we as- 
sume that the dimension of the space-time is D = 4 — 2 e, 
and that of space d = 3 — 2 e. The parameter e is consid- 
ered as small, but only on the level of matrix elements, where 
an analytic continuation to a noninteger spatial dimension is 
allowed. Let us briefly discuss the extension of the basic for- 
mulas of NRQED to the case of an arbitrary number of di- 
mensions. Some basis of dimensionally regularized NRQED 
in the context of hydrogen Lamb shift has already been for- 
mulated in 01, however our approach presented below differs 
in many details. 

The momentum-space representation of the photon propa- 
gator preserves its form, namely g^/k 2 . The Coulomb inter- 
action is 



V(r) = -Ze 2 



i k-r 



d d k 



(2n) d k 2 



Ze 2 



Airr 



1-2 e 



(4 ^ 



r(l-2e) 



r(i 



„l-2e 



(2.1) 



2 



where the latter representation provides an implicit definition 
of Z z , and we have used the formula for the surface area of a 
e?-dimensional unit sphere 
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The nonrelativistic Hamiltonian of the hydrogenic system is 
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(2.3) 



We now turn to relativistic corrections to the Schrodinger 
Hamiltonian in an arbitrary number of dimensions. These cor- 
rections can be obtained from the Dirac Hamiltonian by the 
Foldy-Wouthuysen transformation. In order to incorporate a 
part of the radiative effects right from the beginning, we use an 
effective Dirac Hamiltonian modified by the electromagnetic 
form factors Fx and F 2 (see, e.g., Chap. 7 of |6]), 



p-eFx(V 2 )A + (3m + eFx{V 2 )A 
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where 



3*3 = VM'-V'A*. 
V* = V i = d/dx i ) 



(2.4) 
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Formulas for the electromagnetic form factors Fx,2 can be 
found in Appendix |X| Having the Foldy-Wouthuysen trans- 
formation defined by the operator S (see Ref. |7]), 
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where k = F 2 (0), the new Hamiltonian is obtained via 
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and takes the form 
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The ellipsis denotes the omitted higher-order terms. We adopt 
the following conventions: {X, Y} = XY + YX,Tr=p — 



eA, o % i = [a 1 , a^]/(2i), and the form factors Fx, F 2 are 
defined in Eq. iAl\ below. In d = 3 spatial dimensions, the 
matrices are equal to e y fc a k . The electromagnetic field in 
Hpw is me sum of the external Coulomb field and a slowly 
varying field of the radiation. 

There is an additional correction that cannot be accounted 
for by the Fx and F 2 form factors. It is represented by an 
effective local operator that is quadratic in the field strengths. 
This operator is derived separately by evaluating a low-energy 
limit of the electron scattering amplitude off the Coulomb 
field. An outline of this calculation are presented in Ap- 
pendix[E] The result is 

5H=^E 2 X , (2.10) 
m J 

where E is an electric field, and the functions \ = %W + x' 2 ^ 
are given by Eq. JB12i . 



III. ONE-LOOP ELECTRON SELF-ENERGY 

A. Brief outline of the calculation 

The one-loop electron self-energy contribution in hydro- 
genlike atoms is 
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(3.1) 



Here, p° = E^ is the Dirac energy of the reference state, V 
is the Coulomb potential in d dimensions, and we use natural 
relativistic units with % = c = eo = 1, so that e 2 = Ana. The 
electron mass is denoted by m, and 5m is the one-loop mass 
counter term. By ip we denote the Dirac wave function. There 
are three energy scales in Eq. ( 13. U . which imply a natural 
separation of the one-loop 5^'E into three parts, 



8 {1) E = E L + E_ 
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Each part is regularized separately using the same dimen- 
sional regularization. El is the low energy part, where the 
photon momentum is of order k ~ (Z a) 2 m. Em is the 
middle-energy part, where k ~ m, and the electron momen- 
tum is p ~ [Z a) m. Finally, Ejj is a high-energy part where 
all loop momenta are of the order of the electron mass. It is 
given by the forward three-Coulomb scattering amplitude and 
is represented as a local interaction, proportional to 8 d (r). 

The naming convention for the high-, middle-, and low- 
energy parts is a little different from our previous convention. 
E.g., in |8f], the contribution referred to as the "high-energy 
part" in this reference would correspond to the sum of the 
"high-energy part" and the "middle-energy part" in the con- 
text of the current evaluation. The renaming of the contribu- 
tions is influenced by the NRQED-approach used here and by 
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the correspondence of the different parts to specific effective 
operators. In this work, for all operators Q, we consider only 
the expectation values for states with 



1^0, 



(3.3a) 



and the normalized difference of expectation values 

<(Q)) = n 3 (nS\Q\nS) - (1S\Q\1S) (3.3b) 

for S states. For this reason the high-energy part Eh vanishes 
here. Consequently the "middle-energy part" as considered 
in the current investigation corresponds exactly to the "high- 
energy part" of Refs. 

The one-loop bound-state self-energy, for the states under 
consideration can be written as 

S(1)E= a (Za£ ^ + 2 ^ ln[(Za) - 2] + | 

(3.4) 

where the indices of the coefficients indicate the power of Z a 
and the power of the logarithm, respectively. The coefficient 
A40 is well known (for reviews see e.g. II llll 211 ). and we focus 
here on derivation of the general expression for the a (Z a) 6 
term. 



B. Low-energy part 

In the low energy part, all electron momenta are of the or- 
der of Z a, so in principle, one could perform a direct nonrel- 
ativistic expansion of the matrix element 
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that enters into Eq. ( 13. 11 1. It is more convenient however, in- 
stead of using Eq. i3.l\ . to take the Dirac Hamiltonian with 
an electromagnetic field and to perform this expansion by ap- 
plying the Foldy-Wouthuysen transformation. The resulting 
Hamiltonian, in d dimensions, is given in Eq. (12.9b . Here, we 
can neglect form factors and Hfw becomes (from now on we 
will set the electron mass m equal to unity) 



Hfw 
1 



4 



V • E + a lj (E i t> + Ti j E l ) 



(3.6) 



The contribution from the Coulomb potential V is explicitly 
separated from the additional electromagnetic fields E and B. 
The Hamiltonian in Eq. ( 13. 6> may be used to derive the low- 
energy part which receives a natural interpretation as the sum 
of various relativistic corrections to the Bethe logarithm. We 
use the Coulomb gauge for the photon propagator, and only 
the transverse part will contribute. This treatment of the low- 
energy part is similar to previous calculations 0»l2l]> me dif- 
ference lies in the presence of dimensional regularization. 



The leading nonrelativistic (dipole) low-energy contribu- 
tion is 
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(2n) d 2k 



6 ij 



x(cj> 



E-H-k 



k l w 



<f>) , 



(3.7) 



where by H we denote the nonrelativistic Hamiltonian in d 
dimensions, Eq. ( 12. 3> . The wave function <j>, in contrast to ip 
[see Eq. ( 13. 1H . denotes the nonrelativistic Schrodinger-Pauli 
wave function. In the following, we will denote the expecta- 
tion value of an arbitrary operator Q, evaluated with the non- 
relativistic Schrodinger-Pauli wave function, by the shorthand 
notation (Q). 

After the rf-dimensional integration with respect to k, and 
the expansion in e, Elo becomes |5J] 



E L0 = ( 47 r) e r(l + £ )^ 

SIT 



(3.8) 



'{rUI -Ei{± + ^-1n[2(H-E)]}f 



where we ignore terms of order s and higher. Because the 
factor (4 7r) e r(l + e) appears in all the terms, we will drop 
it out consistently in the low-, middle- and high-energy parts, 
and as well as in the form factors. Moreover, in the two-loop 
calculations discussed below, we will drop the square of this 
factor. The contribution Elq can be rewritten as 

Elo - ^-Za^- £ + b -+H(Za)-^{5 d {r)) 



Aa{Zaf 

3— lnfco- 

Sir n 6 



(3.9) 



where the second term in this equation involves the Bethe log- 
arithm In fco defined as 
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(3.10) 

We consider now all possible relativistic corrections to 
Eq. ( 13.91 . and introduce the notation 



E-H 



. p j 



1 



E-H 



2 (Q 



E-H-k 
1 , 1 



(Q - (Q)) 



(E-H)' P E-H 



(3.11) 



where Q is an arbitrary operator. Sq involves the first-order 
perturbations to the Hamiltonian, to the energy, and to the 
wave function. The first correction E^\ is the modification 
of Elq by the relativistic correction to the Hamiltonian, 



Hr 



+ -Za5 d {r) + -a lJ V l Vp> , (3.12) 
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where S d (r) is a d-dimensional Dirac delta function. One 
could obtain Eli by including this Hr in Eq. A3.9I . How- 
ever, for the comparison with former calculations and for con- 
venience we will return to Eq. d3.7i . and split Eli by intro- 
ducing an intermediate cutoff A 



Eli 




k l k 3 



(3.13) 



We recall the relation A = A (Z a) 2 . In all integrals with an 
upper limit A to be discussed in the following, the divergent 
terms in A will be subtracted. Following earlier treatments 
(e.g., 1 13]), we subtract exactly the term proportional to ln(A), 
but not ln(2 A). The presence of the factor \ under the loga- 
rithm in Eq. (13 . 1 51 is a consequence of this subtraction. 

The quantity 0i can only be calculated numerically. In con- 
stitutes one of three contributions to relativistic Bethe loga- 
rithm C, being defined as in 11311 . 



After the Z a expansion with A = A {Z a) 2 , one goes subse- 
quently to the limits e — » and A — > oo. Under the assump- 
tions ( 13. 3L we may perform an expansion in 1/k in the second 
part and obtain 

£ " = f? Ernr=k*) (3 - 14) 
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(3.17) 



Two others 02, 03 are defined in Eqs. J3.20l > and J3.25i below. 
In this sense, the definition of 0\ in Eq. ( I3.16t corresponds to 
the definition of the low energy part C in Eq. (9) of Ref. 11311 . 

The second relativistic correction El2 is the nonrelativistic 
quadrupole contribution in the conventions adopted in |8j, |9[ . 
Specifically, it is the quadratic (in k) term from the expansion 
of exp(ifc • r), 



After performing the fc-integration and with the help of com- 
mutator relations it reads 
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Here, 0\ is a dimensionless quantity, defined as a finite part 
of the fc-integral with divergent terms proportional to A™ (n = 
1,2,...) and ln(A) dropped out in the limit of large A, 
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In a similar way as for Eli, we split the integration into 
two parts, by introducing a cutoff A. In the first part, with 
the fc-integral from to A, one can set d = 3 and extract 
the logarithmic divergence. In the second part, with the fc- 
integral from A to oo, we perform a 1/k expansion and em- 
ploy commutator relations, with the intent of moving the op- 
erator H — E to the far left or right where it vanishes when 
acting on the Schrodinger-Pauli wave function. In this way 
we obtain 
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Here, 02 is defined as the finite part of the integral [see the discussion following Eq. ( 13. 161 1 
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The third contribution El3 originates from the relativistic corrections to the coupling of the electron to the electromagnetic 
field. These corrections can be obtained from the Hamiltonian in Eq. j2.1> . and they have the form of a correction to the current 



Sf = --p l p 2 + - a 13 k j k-f+- a i] kp 3 - - a ij \7 j V . 
2 2 4 4 



The corresponding correction El3 is 
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We now perform an angular averaging of the matrix element, replace k in the numerator by E — H , and use commutator relations 
to bring the correction El3 into the form 
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We again split this integral into two parts. In the first part k < A, one can approach the limit d = 3, and in the second part k > A 
one performs a 1/ fc-expansion and obtains 
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where (3$ is the finite part of the integral 
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This completes the treatment of the low energy part, which is 

El = Eli + E L 2 + E L 3 



E-H-k' 
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C. Middle-energy part 



We here consider the middle-energy part Em as the contribution originating from photon momentum of the order of the 
electron mass and electron momenta of order Z a. In this momentum region, radiative corrections can be effectively represented 
by electron form factors and higher-order structure functions. Electron form factors F\ and F2 modify the coupling of the Dirac 
electron to the electromagnetic field and the resulting effective Hamiltonian is given in Eq. (12.41 . Here we assume that A ~ 0, A 
represents a static Coulomb potential, and E — —VA° is the electric field of the nucleus. One finds a nonrelativistic expansion 
by the Foldy-Wouthuysen transformation in Eq. d2.8i . and the resulting Hamiltonian [see Eq. ( !2.9bH after putting A = and 
neglecting ^(O) 2 is 



H FW = ^ + eF 1 (V 2 )A°-^-|[F 1 (V 2 ) + 2F 2 (V 2 )](v.S + 2^"^^) 
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The leading {a /it) (Z a) 4 one-loop correction reads 



E M o = (6 {1) V) 



where the "radiative potential" SV is defined as 



SV = 



F{(0) + -F 2 (0) 



(3.27) 



(3.28) 



(3.29) 



and the superscript (1) in Eq. d3.28l > denotes the one-loop component of SV. The expansion of Fi in powers of q 2 is obtained in 
Eq. ( IA3I . Using these results, one obtains 
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Together with the low-energy part Elq in Eq. ( I3.9l l. this gives 
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which is the well-known leading (a/ft) [Z a) 4 contribution to the hydrogen Lamb shift. 

Let us now consider the one-loop correction of relative-order (Z a) 2 . The first contribution Emi comes from the one-loop 
form factors F\ and F2 in Eq. ( I2.9bl > combined with the relativistic correction to the wave function: 



Emi = 2 



(( 




F l[l 


»(0) 




(0) 


16 



V 2 V- 



F 2 (1) (0) 



(T ij VVjT 3 ' 



(E - H)> 



<J|p 2 , V 2 V + 2 tr" V*Vy}\ 



f 2 (1) (Q) 

8 



<(W) 2 ). 



By the superscript (1), we denote the one-loop component of the form factors, as given in Eq. (IA2> in AppendixlAl 
The second contribution Em2 comes from an additional term S^H in the NRQED Hamiltonian, see Eq. ( 12. 101 . 



The corresponding correction to the energy is 



and the total Em contribution is 
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(3.34) 



(3.35) 



D. General one-loop result 



We may now present the complete one-loop correction S^'E up to the order (a/71") (Z a) 6 . It is a sum of the low-energy term 
El given in Eq. P.26> . the middle-energy term Em in Eq. P.35> . and the lower-order term En as defined in Eq. J3.3U . 
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The first two terms corresponds to the a (Z a) 4 term in 
Eq. (13. 4> . whereas the latter terms give the a (Z a) 6 con- 
tribution. The relativistic Bethe Logarithm C, defined in 
Eq. ( I3.17> . consists of a sum of j3x defined in Eq. ( I3.16> . 02 
in Eq. (I3.20> . and 03 in Eq. ( I3.25> . For the convenience of 
the reader we briefly recall that all matrix elements should be 



evaluated in d = 3 spacetime dimensions, which implies 

VV -^4TrZaS 3 {r), (3.37a) 
a-L 

(3.37b) 
(3.37c) 

^VT-tZa^. (3.37d) 
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This concludes the calculation of the one-loop electron self- 
energy. The matrix elements entering into J3.36I are evalu- 
ated below in Sees. ImEllmFl and HTTGl for a number of hy- 
drogenic states and compared to results previously obtained in 
the literature. 



E. Results for D states 

Our aim is to give a few numerical results for some phe- 
nomenologically important hydrogenic states, based on the 
general result ( I3.36l l. For D states, the wave function behaves 
at the origin as ~ r 2 . This means that a few matrix elements, 
such as (V 4 V), are actually vanishing. The following is a list 
of the nonvanishing matrix elements for I = 2: 



nD 



Z a 



1 



(E-H) 



~,P 



nD 



nD 



(Zaf 



Z a 



1118 



<s(> 



7875 25 n 105 n 2 
1 Z a 



(E - H)> 



nD 



nD 



nD 



(Zaf 



(Za) 



1 



1 



_ 2 Z a 



nD 



nD ) = 



709 t 

94500 ~ 150 n + 105 n 2 

(Za) 6 2(n 2 -2) 
n 3 105 n 2 

4 1 
105 ~~ Yr? 



(Zaf 



(3.38a) 

(3.38b) 
(3.38c) 
(3.38d) 



All of the above are evaluated on the nonrelativistic 
Schrodinger wave function. They are finite so that one may 
set the space dimension equal to three. The final results for 
the different fine-structure sublevels are 

A 60 (nD 3/2 ) + A 61 (nD 3/2 ) In [(Za)- 2 ] 

r< n ^ 157 3 , 3007 
= £(nD 3/2 ) 



30240 80?i 37800n 2 



(3.39a) 



and 



A 60 (nD 5/2 ) + A 61 (nD 5/2 ) In [(Za) 



= C(nD 5/2 ) 



379 



1 



1759 



18900 60n 18900n 2 



3l5 ( 3 - 39b > 



They are in agreement with results reported previously in 
Eqs. (12c) and (12d) of 11311 . Values for C(nD 3 / 2 ) and 
C(nD 5 / 2 ) can be found in Table 1 of Ref. II 311 . 



F. Results for P states 

For P states, a few more of the matrix elements in 
Eq. J3.36b are nonvanishing, and we have 



nP 
nP 

nP 

nP 
nP 



Z a 



1 



(E-H) 



IP 



nP 



(Zaf ( 346 4 22 
n 3 \ 135 3n hn 2 
Z a 1 Z a 

nP ) = 



(E-H)' r 3 



{Zaf 



(Za) 



Z a 



nP 
nP 



227 1 1 

540 6 n hn 2 

(Zaf 2 (in 2 -2) 



n 3 15 n 2 

(Zaf (A 

5 15 n 2 
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V 2 [4 7r(Za),5 3 (r)] 
_ (Zaf 8 



nP 



3 ^ n 2 

x3 



(3.40a) 

(3.40b) 
(3.40c) 
(3.40d) 

(3.40e) 



(nPj | i a tj p l [4 7T (Z a) S 3 (r)} p> \ nPj) 

4 (Zaf (l-n 2 ) 



= (nP 



a-L 



nP 



3 n 

The results for the different fine-structure sublevels are 

A 60 (P 1/2 )+A 61 (P 1/2 )\n [(Za)- 2 ] 

r( P ^ 637 1 767 

= C(nP 1/2 ) 



(3.40f) 
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1800 4n 5400 n 2 
29 



15 45 n 2 



\n[\(Zaf 



(3.41a) 



and 



AwCifya) + A 61 (P 3/2 ) In [(Z a)~ 2 ] 

-C( n p ^ J_ 2147 
(n 3/2) + 72()0 + lgn 5400n 2 



(3.41b) 



As for D states, the values for L(nP\j 2 ) and L(nP 3 / 2 ) can 
then be found in Table 1 of fTjll . and the polynomials in n^ 1 
which are part of the above results are consistent with those 
reported in Eqs. (12a) and (12b) of Ref. Ilal . 



G. Results for the normalized difference of S states 

Considering the following matrix elements for / = of the 
S'-state normalized difference ((■)), as defined in Eq. J3.3I >. we 



8 



obtain 

4 7r(Za) S 3 (r) 
= 32 {Zaj 



(E-H) 



7/ 



1 1 5 

4 n 4 n 



+ 7-3 + 7 + *(n) - Inn 



7 4:TT(Za)S 6 (r) 



= 16 (Z a) 



V^[47r(Za) 5 6 (r)}} =8(Za) 



1 h 7 + ^{n) - km 

n 



1-n 2 



(3.42a) 



(3.42b) 



(3.42c) 



(Zaf 



8(Za) 6 



2 1 1 

3 2n 6n 



+ — + 7-1+7 + *(»)- In" 



(3.42d) 



Here, 7 = 0.577216 ... is Euler's constant. One finally ob- 
tains the following result for the general normalized difference 
of the self-energy for S states, 



A 6Q (nS) - A 60 (1S) + [A 61 (nS) - A 61 (1S)} In [(Za)- 2 ] 

rr e ^ rc\Q \ 16087 263 7583 
= C(nS 1/2 ) - £{lS 1/2 ) 



5400 60 n 5400 n 2 



163 



+^[7 + *(n)-lnn]+ln[|(Za)- 



103 4 77 



4[ 7 + *(n) - Inn] } .(3.43) 



45 n A5n 2 



Here, ^(x) = T'(x)/T(x) is the logarithmic derivative of the 
Euler Gamma function. Values for Aqo(tiSi/2) in the range 
n = 1 , . . . , 8 have been obtained using the above formula 
( I3.43> and a generalization of methods used previously for 
states with nonvanishing angular momentum quantum num- 
bers (see Tabled. 



TABLE I: Detailed breakdown of the contributions to A(io{nS), obtained with the help of Eq. J3.43I . The results for Aao{^S) 
and ^60(25) are obtained here with an increased accuracy as compared to Ref. 0]. The generalized Bethe logarithms /3i, P2 
and /33 are defined in Eqs. 13. 1 6i . 13.201 and 13.251 . respectively. The contribution Tt is a contribution to A§o from high-energy 
virtual photons, given in Eq. (5.1 16) of Ref. 1 8] for the 15 state and generalized to arbitrarily high principal quantum numbers 
using Eq. <3~43l We have A 60 (nS) = C(nS) + H(nS) and recall that C = J^Li ft- 



n 


pi{nS) 






£{nS) 


H(nS) 


A 60 (nS) 


1 


-3.268 213 21(1) 


-40.647 026 69(1) 


16.655 33043(1) 


-27.259 909 48(1) 


-3.664239 98 


-30.924 149 46(1) 


2 


-6.057 407 04(1) 


-39.829 658 28(1) 


17.536099 97(1) 


-28.350 965 35(1) 


-3.489499 74 


-31.840465 09(1) 


3 


-6.213 948(1) 


-39.669 430(1) 


17.656 995(1) 


-28.226 383(1) 


-3.476117 


-31.702 501(1) 


4 


-6.167 093(1) 


-39.611903(1) 


17.695 346(1) 


-28.083 650(1) 


-3.478 272 


-31.561922(1) 


5 


-6.100 341(1) 


-39.584 944(1) 


17.712 334(1) 


-27.972 951(1) 


-3.482442 


-31.455 393(1) 


6 


-6.039 851(1) 


-39.570 199(1) 


17.721 349(1) 


-27.888 701(1) 


-3.486429 


-31.375 130(1) 


7 


-5.988 793(1) 


-39.561 272(1) 


17.726711(1) 


-27.823 354(1) 


-3.489 870 


-31.313 224(1) 


8 


-5.946 180(1) 


-39.555 462(1) 


17.730161(1) 


-27.771481(1) 


-3.492 776 


-31.264 257(1) 



One observes the somewhat irregular behavior of (5\ as a func- 
tion of n, which is partially compensated by the other contri- 
butions to Ago. Compared to other families of states with the 
same angular momenta but varying principal quantum num- 
ber 11311 . the A eo for S states display a rather unusual behav- 
ior as a function of n, with a minimum between n = 2 and 
n = 3. The calculations of the relativistic Bethe logarithms 
C, for higher excited S states, are quite involved and will be 
described in detail elsewhere. The value for 15 as reported 
in Table U] represents an improved result (with a numerically 
small correction) as compared to the result communicated in 
Ref. 0], as already detailed in 0- F° r n ^ 3, the results 
for A 60 have not appeared in the literature to the best of our 
knowledge. The results for n = 3 and n = 4 are consistent 



with numerical results for the self-energy remainder function 
as reported in Ref. lTl5ll for these states. 
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IV. TWO-LOOP ELECTRON SELF-ENERGY 
A. Calculation 

The two-loop bound-state energy shift, for the states under 
investigation here, can be written as 



S^E = ^ 

V 7T 



(Za) 4 



{B 



40 



+ {Zaf [B 62 \n 2 [(Z a)' 2 ] + B 61 ln[(Za)~ 2 ] + B 6Q ] } . 

(4.1) 

Here, the indices of the coefficients indicate the power of Z a 
and the power of the logarithm, respectively. The coefficient 
£?4 is well known (for reviews see e.g. fllUl2ll '). and we focus 
here on general expressions for the a 2 (Z a) 6 coefficient. We 
split the calculation into four parts, labeled i — iv according to 
the subsets of diagrams in Figs.^— @] This entails a separation 
of the two-loop energy shift according to 



5 (2) E = s (2) E i + S (2) E U + 5 (2)_g<« + S (2)^ 



(4.2) 



The specific contributions will be considered subsequently in 
the following sections of this article. The B-coefficients cor- 
responding to the subsets i — iv will be distinguished using 
appropriate superscripts. 

We first focus on the pure two-loop self-energy diagrams 
as shown in Fig.^ an d denote the corresponding energy shifts 
and ^-coefficients by a superscript i. As compared to the 
one-loop case treated in Sec. |I11J the two-loop calculation in- 
volves a few more terms with regard to the form-factor con- 
tributions. However, as it has been stressed in Refs. Ildll7ll . 
the leading order of the two-loop low-energy part is already 
(a/n) 2 (Za) 6 , so there are no relativistic or quadrupole cor- 
rections to include at this energy scale. More precisely, we 



split the two-loop contribution into four parts |4|: 

S^E l = E L + E AI +E F + E H . (4.3) 

Here, the contributions El, Em and Ejj are appropriately 
redefined for the two-loop problem [cf. Eq. d3.2l > for the one- 
loop case]. We use definitions local to the current Section for 
the specific contributions. 

The two-loop Ejj is a high-energy part given by a two- 
loop forward scattering amplitude with three Coulomb ver- 
tices. Because it leads to a local potential (proportional to a 
Dirac S in coordinate space), the term Ejj does not contribute 
to the energy of states with I ^ or to the normalized dif- 
ference of S'-states. So, we will not consider this contribution 
here. For S states, this term gives an n-independent contribu- 
tion to the nonlogarithmic term B 60 . 



FIG. 1 : Pure two-loop self-energy diagrams (subset i of the two-loop 
diagrams). The double line denotes the bound-electron propagator. 

The form-factor contribution Ep corresponds to an integra- 
tion region where both photon momenta are of the order of the 
electron mass, but the electron momentum is of the order of 
Z a. This part is a sum of two terms: 



Ep — Epi + Ep2 ■ 



(4.4) 



The first term Ep\ comes from two-loop form factors, in the 
same way as the one-loop Emi [see Eq. J3.32H . It contains 
additionally an iteration of the one-loop potential S^'V and 
the term proportional to k 2 from Eq. ( I2.9bl i: 



E F1 = (S^V*) + 2 (s^V* 1 H R ^ + f i (2) (°) + 2J1 2 (2) (0) { ^4y + 2 ia i *p i V 2 Vpl) + F^ 2 \o) (V 4 V) 



F 2 (2) (0) 



1G 



{p 2 , V 2 V + 2a lJ V l Vp 3 Y) 



F 2 (2) (0)+ Lf 2 (1) (0) 



8 



((w) 2 > + <Vv^-^mv 



(4.5) 



The two-loop form factors are given in Eq. l|A3i below, and S^V, 5 {2) V are the one- and two-loop components respectively of 
the potential given in Eq. J3.29I . The explicit form of 8^ 2 ' V % can be found in Eq. j4.20t below. 

Ep 2 comes from the low-energy two-loop scattering amplitude and is the analog of the one-loop Em2 m Eq. ( I3.34> . The 
effective interaction is 

S (2) H = x (2) ^ V) 2 > (4 6) 

where is defined in Eq. JB12bi below. It is assumed that vacuum polarization diagrams does not contribute in the current 
section to form factors as well as to x- The energy shift due to S^'H is 

E F2 = (S^H\ . (4.7) 

It is a remarkable fact that this two-loop scattering-amplitude contribution is infrared finite, in contrast to the corresponding 
one-loop result in Eq. (I3.34> . 
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For the two-loop problem, we redefine Em to be the contribution where one of the photon momenta is of the order of the 
electron mass, the second photon momentum is of order (Z a) 2 and the electron momenta are of order Z a. In the spirit of 
NRQED, the contribution coming from large photon momenta is accounted for by form factors. Therefore Em is given by the 
correction to Bethe logarithms coming from one-loop form factors. It is a sum of two parts 



Em — Emi + Em2 ■ 

The contribution Emi is similar to the one-loop term Eli with Hr replaced by S^V: 



E 



Ml 



d d k d-1 



%i>V ( P 



1 



■P 



{2n) d 2k d "°"" v Y E-H-k 
We calculate it by splitting the integral in two parts k < A and k > A in analogy to the one-loop case, 

1 



(4.8) 



(4.9) 



E 



Ml 



., f A d d k d-1 . 

■ dsn 



E-H-k 



r d° 

" J {2n) d 2k ~d 

+ ^{( K[ ,av, p1] ) + 2 (,aV^v 



v 



(4.10) 



where 



= h + { l~ \ ln[2 {z a)2] } + £ { I \ c(2) + t ln [h {z ay2] + \ w [ ^ {z ay2] } ■ (4 - 1 1} 

We have not approached the limit d — 3 in the first part, because S^V contains 1/e. It will eventually cancel when combined 
with El, and only then one approaches this limit. Emi is similar to the one-loop Els and comes from the ^(O^correction to 
the coupling with the radiation field, 



H FW = -- a 11 [1 + F 2 (0)] - | [1 + 2 F 2 (0)] [V • E + ^ (E l tt j + ir j E l 
4 8 



which yields 



W = 2 o a ' k (k-rk k +ikp k - V k V) 



The corresponding correction Emi is 



E 



M2 



2e z 



d d k 



(2ir) d 2k 



gij 



k l W 



a 2d-3 
"n 4 (d-1) 



5f 



E-H-k 



a ik yky 



(4.12) 



(4.13) 



1 fa\ 2 
2 



dkk ( a %d VW 



E-H-k 



P 



(4.14) 



and this integral in (Z a) 6 order does not depend on the cut-off in the limit A — > oo, when one drops the linear term in A. 

The low-energy part El, appropriately redefined for the two-loop problem, is a contribution from two low-energy photon 
momenta, fc; ~ (Z a) 2 . Its explicit expression is rather long: 



Er 



I — 

J (27T 



d d ki d - 1 



P(k 1 ,k 2 ) = (p i 



) d 2h d 
1 



d d k 2 d - 1 



(2n) d 2k 2 d 
1 



P(k u k 2 ), 
1 



■ p J 



■P 



E-(H + h)" E-{H + ki + k 2 ) r E-(H + k 2 ) 
1,1 1 



+ 2 \ P E- (H + k 1 ) t ' E-(H + k 1 +k 2 ) r E-(H + h 



■p> 



E-(H + k 2 )^ E- {H + ki +k 2 ) E-(H + k 2 ) 



+ {P 



E-{H + k 1 ) t ' (E-H)~ P E-(H + k 2 ) 



1 



P 



P ) (P 



1 



■ p J 



1 



E-(H + k 1 ) 1 ' / Y [E-(H + k 2 )} 2 
+ (pS fir , , ^ 77T-, , \ Pf - , 1 , \ P 



1 



P 



2 y E-(H + k 2 

1 ,A I 

P 



P ) (P 



1 



■ p 3 



E - (H + ki) E - (H + k 2 ) I h+k 2 Y E-{H + k 2 y / k 1 + k 2 \ E - (H + fci) 



[E-iH + h)} 2 

p- ,1 , , y )- (4.i5) 
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We calculate El by splitting both integrals in a way similar to the derivation presented in |4], 



E, 



a\ 2 (Zaf , 2 
1 b L + e 2 



d d k 2 d - 1 a £ 



(2 7r) d 2A:2 d tt 2 



<$v 2 v ( P 



E-H-k' 



-a f 


2 r 


_7T 2 





/vV 

\ {E-H)> 



vV 



Here, the two-loop Bethe logarithm bp is obtained as the finite part of the integral 



(Zaf 4 



n" 



dk\k\ I dk 2 k 2 P(ki, k 2 ) . 



where it is assumed that the following limits are performed in order: first d — > 3, next A2 — > 00 and finally Ai 
above. This definition of bp corresponds to the one in Refs. II all 711. 



(4.16) 

(4.17) 
00 in the 



B. General result for the pure two-loop self-energy 



The pure two-loop self-energy contribution up to the order a 2 (Z af, denoted 8^E l (see Fig.Q, may now be obtained as 
the sum of Ep + Em + El- With the partial results given in Eqs. ( I4.4L ( 14. 8> and ( 14. 151 . respectively, we obtain 



5 (2) E i = (gWyi^ + 



a\ 2 (Z a) 



6 r 



+ 2 ( 8 {2) V l 



1 



-H R 





10 

¥ + 




' 31 


3 




256 


+ 16 



Vj 



(E - H)> 



Vi 



a N 

7T- 

a\ 2 

71". 



(E-H) 

S + yC(2)M2)-|«2)-g«3) 



5 

32 



3 

64 



{i, 



}) 



<(W) S 



3295 9 



4063 



777^) ln(2i vr^C(2) - ^C(3) + |r In [i (Za)- 2 } + 1 In 2 [± (Za)" 2 ] 



41472 10 

3059 1 , 

- - C(2) ln(2) 



14400 
1321 



54 



18 



(v¥> 



23040 5 

Here, the first term (S^'V 1 ) is of lower-order [a 2 (Za) 4 ], and 

V 2 V" /10 4 



f(2) + — C(3) + — In [i (Z a)' 2 ] 
5760 sv ; 20 sw 24 L2 1 ; J 



(2ia ij p* V 2 Vp>). 



(4.18) 



Vi 



a 

7T 

(- 

V 7T 



4 V T + 3 ln ^ (Zar2 ] + X v ^ pJ 



163 9 , „ 85 9 ,.„ x 

+ - C(2) In 2 C(2 C(3) 

288 4 ^ W 144 SV ; 16 ^ W 



(4.19) 



V 2 V 



-|-|«2)ln2 + 5c(2) + 2c(3) 



(4.20) 



The various generalized Bethe logarithms that enter into Eq. J4.18L are given as follows (with the implicit assumption that 
polynomial divergences as well as logarithmic ones for large A = A/(Z a) 2 are dropped) 



(Zaf 



(Zaf 



(Zaf 



N = —Za / dkkS v g3/ r \(p 
3 Jo 



/?4 = — / dk kSt 



1 



3 Jo 



(o i iV i Vpi/i) \ P 



E-H-k' 
1 



E-H-k 
1 



^ = U dkk {-r l ' V3V E-H-k 



p 



(4.21a) 
(4.21b) 
(4.21c) 



The N term has previously been defined in Refs. 141 11611 : it is generated by a Dirac delta correction to the Bethe logarithm. All 
the explicit matrix element occurring in the formula (14. 1 8I > can be calculated using standard techniques, for arbitrary hydrogenic 
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states with nonvanishing angular momentum, and for the normalized difference J3.3I > of S states. The evaluation of the general- 
ized Bethe logarithms N, P4, and (3$ is more complicated (see Refs. 1 13l fl8l1 ~). The calculation of the two-loop Bethe logarithm 
£>l for arbitrary excited hydrogenic states is a challenging numerical problem. So far results have been obtained only for excited 
S states U7I1 . The formula (14. 1 8i thus provides the basis for complete two-loop calculations in the order a 2 (Z a) 6 , and 
reduces the remaining part of the problem, for a general hydrogenic state, to a well-defined and in essence merely technical nu- 
merical calculations. In the following sections we discuss the evaluation of the formula (I4.18> for particular hydrogenic states for 
which the generalized Bethe logarithms can be inferred from previous calculations. These comprise the fine-structure difference 
of D and P states, and the normalized difference for S states. 

TABLE II: Numerical values for the pure two-loop self-energy diagrams as shown in Fig.Q The Beo -coefficients receive a superscript i. 



n Bi (D5/2-D 3/2 ) BUP3/2 - P1/2) b L {nS) N(nS) R{n) B^nS) - B^IS) 



1 — — -81.4(3) 17.855 672 03(1) — — 

2 — -0.361 196 -66.6(3) 12.03214158(1) -0.671347 14.1(4) 

3 -0.018 955 -0.410 149 -63.5(6) 10.449 809(1) -1.041532 16.9(7) 

4 -0.022 253 -0.419 927 -61.8(8) 9.722 413(1) -1.254 980 18.3(10) 

5 -0.023 395 -0.420 828 -60.6(8) 9.304114(1) -1.392 573 19.4(11) 
6 -0.023 826 -0.419 339 -59.8(8) 9.031832(1) -1.488 456 20.1(11) 



C. Results for the fine-structure difference of D states 

For D states, we use the general result (14.181 and the fact that matrix elements involving a Dirac-delta function vanish. Thus, 
logarithmic terms for D levels vanish, B^^D^m) = Bq 1 (D 5 / 2 ) — 0. The absence of logarithmic terms even holds for the sum 
of all two-loop diagrams (not only for the subset i), and even for arbitrary states with orbital angular momentum I > 2. This 
result generalizes the well-known fact that the double-logarithmic contribution B§ 2 vanishes for states with / > 2 l!9H2f)ll . For 
the fine-structure difference of _B| , we use the result in Eq. ( 14.181 1 and the matrix elements in Eq. ( 13.381 . to obtain 

38497 133 895 / 3817 13 29 \ 
B 6Q (D 5/2 D m ) - - — + + - — + _ j C (2) ln(2) 

/ 3817 13 145 \ / 3817 13 29 \ , /0 , „ „ . „ „ . 

+ (3024O + 48^ - 168^) C(2) + (j00800 + 160^ - TT2^J C(3) + WW ~ D m) + W*/2 - ■ 

(4.22) 

Numerical data for Bqq(D 5 / 2 — -D3/2) can be found in Tab. ITT1 The unknown two-loop Bethe logarithm 6i(nD) does not 
contribute to the fine-structure difference of D states. 



D. Results for the fine-structure difference of P states 

We again use the fact that the unknown two-loop Bethe logarithm &l does not contribute to the fine-structure difference of P 
states. With the help of the general result in Eq. (14. 1 8i and the matrix elements in Eq. ( I3.40> . we obtain 

S 61 (P 3/2 -P 1/2 ) = ~ (l-^) (4.23) 

in agreement with the literature l2lll and 

, . , 217 151 325 (\ 1 \ , . . / 103 15 37 \ . , . . 

Bl (P 3/2 - P 1/2 ) = -— - — + — 2 + (- - — 2 ) ln(2) + {-— - - + _ j C (2) ln(2) 

23 25 749 \ , N /103 15 37 \ , N 

160 + T6^" 720^) C(2) + (960 + 32^ - 80^J C(3) + ^ ~ P ^ + ^ ~ P ^ ' (4 ' 24) 

Numerical values of the relevant quantities for n = 2, ... ,6 can be found in Tab. ITTI They are in full agreement with results 
previously obtained in l2lll . The generalized Bethe logarithms (3^ and f] 5 in these expressions are equivalent to the quantities 
Af s ^4(n) and Af s £ 5 (n) as defined in Ref. 12111 . In the context of the current investigation, the numerical values of Af s £ 4 (n) and 
Af s £s(n) were reevaluated with improved accuracy as compared to Ref. \22\ and the data in Tab.[I]]are consistent with them. 
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E. Results for the normalized difference of S states 

We evaluate the general formula given in Eq. d4 . 1 8 i for the normalized difference of S states, using the matrix elements given 
in Eq. ( 13.42k In the result, we identify terms with the square of the logarithm ln[(Z a)~ 2 } (B\ 2 coefficient), and with single 
logarithm (B l &1 coefficient), and the nonlogarithmic term _B 60 . The results discussed here probably are the phenomenologically 
most important ones reported in this paper, because of the high accuracy of two-photon spectroscopic experiments which involve 
S—S transitions. 

For the double-logarithmic term, we recover the following known result (see Refs. H E^ V 
Bl 2 (nS) - Bi 2 (lS) = ^0 + 2 L-i+ 7 -l n (n)+*(n) 



Here \& denotes the logarithmic derivative of the gamma function, and 7 = 0.577216 . 
restricted to the two-loop diagrams in Fig.Q] reads I4H23I1 



(4.25) 

is Euler's constant. The result for Bqi, 



B l 61 (nS) - B* 61 (1S) = - [N(nS) - N(1S)] + 



80 
27 



32 
~9~ 



In 2 



1 



+ — it + 7 - ln(n) + *(n) 



(4.26) 



This result is recovered here from Eq. ( 14.181 1. using the matrix elements in Eq. ( 13.421 1. Moreover, we obtain the complete 
n-dependence of the nonlogarithmic term £>g : 

BUnS) - B\ Q {IS) = b L {nS) - b h {\S) + - 5 m(2)) WnS) - N(1S)} 



10529 14099 



17699 



5184 
53 
15 



2592 n 5184 n 2 
35 149 
2n 



4 16 4 . , , 
3-9^ + 9^) ln(2) 



20 
~9~ 



80 
27^ 



20 
27^ 



30 rv 



C(2) ln(2) 



497 16 



/1357 167 2792 
V2700 ~ 36^ + 675 n 2 
79 



C ( 2 ) + [ !3 _ 35 
^ K ' 1 60 8n 



ln(2) 
419 



120 rv 



C(3) 



1296 
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The generalized Bethe logarithms j3i and /3s, which make an occurrence in Eq. ( I4.18t but are not present in Eq. (14.271 . vanish 
for S states. The result J4.27i can also be written as 



B l 60 (nS) - B 1 60 (1S) = b L (nS) - b L (XS) + R(n) 



(4.28) 



which provides a definition of the remainder R(n). Numerical values for b^inS), N(nS), R{n) and the normalized S-state 
difference B l &Q (nS) — Bq Q (1S) are given in Table UTI and we have the opportunity to correct a calculational error for N(2S) 
whose value had previously been given as 12.032209 in fUjil . 




FIG. 2: Feynman diagram with a vacuum-polarization loop in the 
self-energy virtual photon line (this single diagram forms subset ii 
in the convention adopted in this paper). 



V. FERMION LOOP IN THE SELF-ENERGY PHOTON LINE 



We here calculate the mixed self-energy vacuum-polarization diagram in Fig.|2] The result can be easily inferred from the 
terms in square brackets in Eqs. (14. 6> . ( IA3I . and JB 12i . and reads 
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where 
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"324 + 144^ ' 
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(5.2) 



is a radiative potential in the sense of Eq. J3.291 . but includes here only the vacuum polarization part of form factors. We observe 
the absence of ln(Z a) terms. 



TABLE III: Values of the fine-structure for D and P states and normalized difference of 5* 
states coming from from fermion loop diagrams in Fig. [2] 



n B&(D t/a - D B/a ) B&(P B/a -P 1/2 ) B&(nS) - B&(1S) 

2 — -0.013 435 0.109 999 

3 0.000 757 -0.017 089 0.114 502 

4 0.000 878 -0.018 613 0.110 743 

5 0.000 915 -0.019 431 0.106 566 

6 0.000 925 -0.019 935 0.102 982 



Numerical values for S, P, and D states can now be obtained using matrix elements in Eqs. (13.381 . (13.401 and J3.42t . For the 
fine-structure intervals, we obtain 



B&(D 6/3 ~ D 3/2 ) 
B eo(P3/2 - Pi/ 2) 
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Considering S states, as is evident from Eq. i5.1\ . using the matrix elements in Eq. 13.421 . the normalized difference of £>g\ 
vanishes, Bq X (nS) — E>1\ (15) = 0, and this result is in agreement with the literature. For the normalized rt-dependence of Bg , 
we obtain the following result, 



B$ s (nS)-B&(lS) = - 
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Numerical values are presented in Tab. Mil 



VI. COMBINED SELF-ENERGY WITH A FERMION LOOP IN THE COULOMB PHOTON LINE 

The Feynman diagrams in Fig. [5] represent the modification of a leading one-loop self-energy correction by a perturbing 
Uehling potential Vu = ~jg a {^ a ) 5 d (r) — V 2 V. One can easily obtain the result from Elq + Emo m Eqs. ( 13.91 and 

( I3.30> . by replacing the Coulomb potential V by V + Vu, and expanding all matrix elements in Vu, up to the linear terms. The 
result is 
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(6.1) 



Here, N is a correction to the Bethe logarithm as defined in 
Eq. ( 14.2 la^ All matrix elements in this result vanish for in- 
states and for states with higher angular momenta. The ab- 
sence of both logarithmic as well as nonlogarithmic terms 



holds from subset Hi holds for arbitrary states with orbital 
angular momentum / > 2. For P states, we obtain the fine- 
structure difference i?g\ l (nP 3 / 2 ) — Bq{ (nPi/ 2 ) = 0, in agree- 
ment with the literature. For the nonlogarithmic term, we ob- 
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tain 



is 



B%j(nP 3/2 ) - B™(nP 1/2 ) = Lfi-^\. (6.2) 

As a last example, we consider the 5-state normalized dif- 
ference defined in Eq. ( 13 .31 . using matrix elements given 
in J3.42i . For the double-logarithmic term, we recover the 
known result BjM(nS) = (see Refs. HHj). The result for 
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Bfl reads HQ 
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where the one- and two-loop parts read 1 24, 25] as follows, 

t 2 \ 4- ,7,( 2 ) (-n 2 \ 4- 



15 140 
41 449 <f 2 
162 ~ 10800 
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(7.3) 



This result is recovered here from Eq. ( 16. 1> . As a new result, 
we obtain the complete n-dependence of the nonlogarithmic 

3 60 ■ 



B$(nS) - B%(lS) = -- [N(nS) - N(1S)} 



16 



9 27n Tin 2 
16 32 
27 + 45 



32 



r)M2) 



15 45n 45n 2 , 
ln(2) ) [7 + *(n) - ln(n)] . (6.4) 



Using this formula, it is then possible to infer the values of 
B$(nS) - B™{1S) as given in Tab.Hvl 



In the integral representation for lo^ given in Eqs. (15) 
and (16) of Ref. |8|. one should make the replacement 

In fj^f J —> m (l^f) m (^r) in order t0 correct f° r a ty- 
pographical error in an intermediate step of this calculation. 
In the coordinate space, the correction becomes 
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The contributions to the energy involves the first and second 
order matrix element together with relativistic corrections, 
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FIG. 3: Two-loop diagrams (subset Hi) generated by a fermion loop 
in the Coulomb exchange of a one-loop self-energy. 



TABLE IV: Values of the difference B%$(nS) - B$(1S) for the 





n 


Bg}(nS)-B$(lS) 


2 


1.491 199 


3 


1.890 577 


4 


2.072 903 


5 


2.177 348 


6 


2.245 177 



VII. PURE TWO-LOOP VACUUM POLARIZATION 

We investigate the subset of Feynman diagrams in Fig. |4] 
The vacuum polarization correction to the Coulomb potential 



The two-loop part of this expression reads 

41 /a\ 2 
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(7.6) 



The first term in this result corresponds to the a 2 (Z a) 4 term 
in Eq. J4-. li . The remaining terms give the B^ coefficient. 

We first notice the complete absence of logarithmic terms 
in the result ( 17. 6> . All matrix elements in (17. 6> vanish for in- 
states and for states with higher angular momenta, in the order 
of a 2 (Zaf . The fine-structure difference of the nonlogarith- 
mic term for P states is as follows, 

B™{nP m ) - B%(nP 1/2 ) = (l - -L\ (7.7) 
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The n-dependence of the nonlogarithmic term B£q is as fol- 
lows, 



B&(nS)-B&(lS) 



1817 2194 1337 
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This completes our investigation of the subset iv. 



(7.8) 



FIG. 4: The remaining two-loop diagrams (subset iv) involve at 
least one closed fermion loop in the Coulomb photon exchange be- 
tween electron and nucleus, and no self-energy photons. 



TABLE V: Values of the difference B x a v {nS) - B^IS) for the 
diagrams in subset iv. 



n Bl v (nS) - Bjg(lS) 



-0.611365 
-0.603 468 
-0.560 004 
-0.521300 
-0.490 240 



VIII. TOTAL RESULT FOR ALL TWO-LOOP DIAGRAMS 



The two-loop subsets i — iv (see Figs.^— E} have been con- 
sidered in Secs. lIVI — IVIII We are now in the position to add 
the results given in Eqs. ( I4.18> . ( 15. lL ( 16. U and il.6t , and to 
present a general expression for the complete two-loop cor- 
rection to the Lamb shift, including the vacuum-polarization 
terms, valid for general hydrogenic bound states with nonvan- 
ishing angular momenta, and for the normalized difference of 
S states. This general result reads 
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The third line in the above equation corresponds to the lower-order a 2 (Za) A contribution (B40 coefficient). We now turn to the 
evaluation of this expression for S states. The sum of the contributions in Eqs. 14.27L J5.5I >. ( 16.41 and fl7.8i corresponds to the 
sum of all the matrix elements in Eq. (18. 1> . evaluated for the normalized difference of S states. The logarithmic terms 0] have 
already been verified for the normalized difference. The n-dependence of the total nonlogarithmic term may be expressed as 

B 6Q {nS) - B 60 (1S) = b L (nS) - b L (lS) + A(n), (8.2) 

where A(n) is an additional contribution beyond the n-dependence of the two-loop Bethe logarithm, defined in analogy to 
Eq. ( 14.281 . The result for A is 

,/ \ /38 4, .A rAr/ ON AT ,, mi 337043 94261 902609 (A 16 4 \ 2 .. 
A{n) = (l5 - 3 ln(2) ) [NinS) - N(1S)] - TmOO - 21600^ + 129600^ + (3 " 9^ + 9^ ) ^ (2) (8 " 3) 

76 76 ^ m(2) + f-S + £ - ^ C(2) ln(2) + (^-^ + ^) f (2) 



45 135 n 135 n 2 / w V 15 2ra 30 n 2 / w V 10800 2n 10800 n 2 

'I - I + 15Sf ) « 3 » + S + T '° 2(2) - H '° (2) + 8 « 2) ln(2) - f « 2) - 2 « 3) ) h + *<"> - ta< "" 



Numerically, A(n) is found to be much smaller than b^nS) — 
b]j(lS), as shown in Table IVll This implies that the numer- 
ically most important contribution to B 60 (nS) — B eo (lS) is 
exclusively due to the two-loop Bethe logarithm. The theo- 
retical uncertainty of BQo(nS) — Bqq(1S), for higher excited 
nS states, is caused entirely by the numerical uncertainty of 
the two-loop Bethe logarithm bL(nS), with explicitly data for 
higher excited states taken from Ref . | jjj • 



TABLE VI: Total values of the difference B 60 (nS) - B 60 (1S) 
coming from all diagrams. 



n 


A(n) 


B 6 o(nS)- J B 60 (LS') 


2 


0.318 486 


15.1(4) 


3 


0.360 079 


18.3(7) 


4 


0.368 661 


20.0(10) 


5 


0.370 042 


21.2(11) 


6 


0.369 462 


22.0(11) 



For the fine-structure difference of D states, the total two- 
loop results is obtained by evaluating the general result in 
Eq. J8 . 1 1 on D states, or alternatively by adding just the con- 
tributions from subsets i and ii [see Eqs. ( I4.22i and J5.3I 1. 
because the subsets Hi and iv do not contribute to the D fine 
structure. For the P-state fine structure, the sum of the results 
in Eqs. ( I4.24L ( 15. 4t . ( 16. 2t , and J7.7I > gives the complete result, 
including the nonlogarithmic term Bqq. It has already been 
stressed that in order to determine the absolute value of Bqq 
for P and D states, an evaluation of the Bethe logarithm 6l for 
these states would be required, and its knowledge is currently 
restricted to S states. 

Despite this, we may evaluate general logarithmic terms for 
P and D states. For D states and states with higher angular 
momenta, a direct evaluation of Eq. J8. 1 1 immediately reveals 
that the logarithmic terms vanish, 

B 62 (nD) = B 61 (nD) =0. (8.4) 



The same holds for any hydrogenic states with orbital angular 
momentum / > 2. For P states, an evaluation of ( 18. It confirm 
that 

4 n 2 — 1 

B 62 (nP) = (8.5) 
27 n z 

Furthermore, the logarithmic terms are 

(8.6) 

4 T? 2 - 1 / 31 8 \ 

B 6 i(nP 3/2 ) =-N(nP) + ^— — - — In 2 . 

by 3/27 3 V ' n 2 V405 27 / 

(8.7) 

Numerical values for N(nP) can be found in Eq. (17) of 
Ref. 0. 

IX. SUMMARY 

We have presented a unified approach to the one- and two- 
loop electron bound-state self-energy correction in hydrogen- 
like atoms, including terms of order a (Z a) 6 and a 2 (Z a) 6 , 
respectively. We consider states with nonvanishing orbital an- 
gular momentum and the normalized difference of S states. 
The general analytic structure of the one- and two-loop cor- 
rections is given in Eqs. J3.4I and ( 14. 1> . respectively. The gen- 
eral result for the one-loop correction is given in Eq. ( 13.361 1. 
We evaluate our formulas for specific families of hydrogenic 
states in Sees. MI El MI Fl and MI Gl (one-loop case). All one- 
loop results are in agreement with those previously reported 
in the literature. In addition, we obtain results for the nonlog- 
arithmic terms (Aqo coefficients), for higher excited S states, 
as listed in Tab.|J 

For clarity, we separate the two-loop calculation into four 
different subsets i, ii, Hi and iv consisting of separately 
gauge-invariant diagrams (see Sees. IIVI — IVIll and Figs. \l\ — 
[4}. A general formula for the "pure" two-loop self-energy 
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diagrams is presented in Eq. (14. 1 8i . The corresponding ex- 
pression, for the self-energy vacuum-polarization diagram in 
Fig. 13 can be found in Eq. (15- II . For the subsets Hi and iv, 
we present general expressions in Eas. l6.1l and 17.61 For the 
total sum of the two-loop effects, a summary is provided in 

Sec.rvTnl 

The two-loop fine-structure difference for P states for the 
subset i as given in Eq. ( I4.24l > is in agreement with previous 
results l2llE3l . This constitutes an important cross-check of 
the method used in the current investigation, which is based 
on dimensional regularization, and on effective operators for 
the contributions stemming from hard virtual photons. The 
results given in Eqs. J5.4I >. \6.2l . and M .11 complete the fine- 
structure difference of P states in the order a 2 (Zaf . 

The central result of the current investigation, however, is 
the complete n-dependence of all two-loop logarithmic and 
nonlogarithmic contributions to the Lamb shift of S states 
up to the order a 2 (Z a) 6 . In this regard, our study follows 
a number of previous investigations on related subjects (see 
Refs. I4l 12a 1271 12810 . where the logarithmic terms were pri- 
marily investigated, but the nonlogarithmic term was left un- 
evaluated. The n-dependence of all logarithmic terms for 
S states [corresponding to the Bq 2 and Bqi coefficients in 
Eq. ( 14. 1H is recovered in full agreement with the literature. 
For the B§i coefficient, we refer to Eqs. \A.2bl and ( I6.3I >. 
Moreover we obtain in Appendix[C]an additional logarithmic 
contribution B§i (15) to the ground IS state, which was omit- 
ted in the former work pj] . 

Partial results for the n-dependence of the nonlogarithmic 
term £>6o ( n S) are given in Eqs. 14.271 . 14.271 . ( 16. 4t and ( I7.8l l. 
A summary including all two-loop subsets is provided in 
Eqs. ( 18. II . J8.2t and ( 18 .31 . Our results lead to predictions for 
the 5-state normalized difference with an accuracy of the or- 
der of 100 Hz (see Ref. jH and Appendix0. We find that the 



largest contribution to the n-dependence of Bqq stems from 
the two-loop Bethe logarithm b^, but the remaining contribu- 
tions in Eqs. ( 15.51 . ( I6.4I > and. ( 17. 8> are essential for obtaining 
complete predictions (see also Tables IVTl and I Vllt . 
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APPENDIX A: ELECTROMAGNETIC FORM FACTORS 

We consider the form factors defined by 

1,^ r M = FxtfYtp + ^F 2 (q 2 ) f~) , (Al) 

where q is the outgoing photon momentum. The form factors 
are expanded in a up to the second order, 

F^q 2 ) = l + ifVHifV), (A2a) 
F 2 (q 2 ) = F^(q 2 )+F^(q 2 ), (A2b) 

where the superscript corresponds to the loop order, i.e. to the 
power of a. They have recently been calculated analytically 
by Bonciani, Mastrolia and Remiddi in [31]. The results for 
the form factors expanded into powers of q 2 up to q 4 read (in 

D = 4 - 2s): 
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The subscript VP denotes the contribution to the two-loop 
form factors which involves a closed fermion loop (see Fig.[2ji- 



We can now turn to the scattering amplitude T. The expres- 
sion corresponding to the tree diagram of Fig.|5]is 



APPENDIX B: LOW-ENERGY LIMIT OF THE 
SCATTERING AMPLITUDE 

In the leading order, the electron self-energy can be incor- 
porated by electromagnetic form factors F\ and F 2 , and more 
precisely by the leading terms of its low momentum expan- 
sion. In the higher order, namely a (Z a) 6 , single vertex form 
factors Fi are not sufficient, and the additional term is the low- 
energy limit of the spin-independent part of the scattering am- 
plitude with two 7 vertices, (see Fig. [5}, with the form factor 
contributions subtracted. This term, has not yet been consid- 
ered in the literature. A detailed derivation is postponed to a 
separate paper; here we present only a brief derivation. 
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► T - t - 



= I Tr 
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70 (#3 + 1) 
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and this expression defines our normalization. The presence 
of 70 in Eq. (IB6> results from the fact, that we consider the 
scattering by the Coulomb potential 



e^A^ = 70 V = -70 



Ze 2 



(B7) 



The momenta pi and p$ are on mass shell (p\ — 1, p 2 = 1). 
Let us define the exchange momenta according to 



qi = pi - P2 



92 = P2 - P3 ■ 



(B8) 



and the static momentum t, such that t — (1,0) and t 2 = 1. 
Because we consider the scattering of a static potential, the 
exchange momenta are spatial, 



FIG. 5: Tree and one-loop diagrams with two Coulomb exchanges. 

To construct the projection operators for a spin independent 
part of the scattering amplitude, let us consider the matrix el- 
ement of an arbitrary operator Q, namely u(p', s') Qu(p, s), 
where u(p, s) is a positive solution of the free Dirac equation, 
normalized according to uu = 1. We transform this matrix 
element to the more convenient form 



i(p',s')Qu(p, s) = Tr[Qu{p, s)u(p',s')} 



(Bl) 



Because we aim to calculate only the low energy limit, we can 
use an approximate form of u(p, s), 



(B2) 




(B3) 



where I is the 2x2 unit matrix, the spin-averaged projection 
operator becomes 



^2u{p,s)u(p',s) 
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The spin-averaged matrix element of an arbitrary operator Q 
can now be expressed as 
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The one- and two-loop radiative corrections, and T^ 2 \ 
are obtained using standard rules of quantum electrodynam- 
ics. However, we additionally subtract from these amplitudes 
the corresponding form factor contribution. This subtraction 
is carried out using the tree diagram with the vertex 70 re- 
placed by To, 



(A + i)r (<h) 



1 
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r (<?2 



1: 
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The vertex function T M is defined in Eq. dAli . In the one-loop 
order, the subtraction permits the approximation T M « 1 for 
one of the vertices, with a form-factor correction at the other, 
and a second term where the approximations at the vertices are 
interchanged. For the two-loop case, it is understood that the 
subtraction includes only (a/ir) 2 terms, so there are a total 
of three terms, one with both vertices modified by one loop 
corrections, and two others where only one vertex receives a 
two-loop correction. After the form-factor subtractions and 
small momenta expansion, the scattering amplitude takes a 
simple form 



9i ■ 92 X 
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where the superscript denotes the loop order. The coefficients 
X have been calculated with the help of the symbolic program 
FORM HI and read 
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where the subscript vp denotes the contribution from the di- 
agram in Fig. |3 Using the relation qi ■ q 2 = —q\ ■ q\, and 
including the factors given by the Coulomb potential, one ob- 
tains the effective interaction Hamiltonian 



5H 



-(-iW) 2 X = {VV) 2 x = e 2 E 2 X . 



(B13) 



where by — > we denote the transition to the coordinate space 
by the corresponding Fourier transform. 



APPENDIX C: ADDITIONAL LOGARITHMIC 
CONTRIBUTION TO THE GROUND STATE LAMB SHIFT 

The two-loop logarithmic contribution to the Lamb shift has 
been considered by one of us (K.P.) in 0]. The obtained re- 
sults for Bqi coefficient of the ground state was 



^61 



39751 4 . m 
10800 + 3 ^ 



55 7T 2 616 1n(2) 



27 



135 



3 7r 2 ln(2) 40 In 2 (2) 9 £(3) 



8 



50.309 654. 



(CI) 



After careful reanalysis of the performed calculations we 
found that there is an additional logarithmic contribution, 
which can be associated to the e 2 E 2 = (WV) 2 term in the 
effective Hamiltonian in Eqs. J2.9bl2.10l 



5H 



(1)\ 2 



/(2) 



(W) 2 



(C2) 



Since this additional contribution is numerically small, it does 
not explain the discrepancy with the direct numerical calcula- 
tion by Yerokhin et al. in Ref. yD, although the difference is 
now slightly smaller. We postpone further conclusions until 
the evaluation of the constant term i? 6 o is completed. 



APPENDIX D: EVALUATION OF THE LAMB-SHIFT 
DIFFERENCE 

We denote the Lamb shift of an nS states by AE(nS) and 
use the definition in Eq. (67) of Ref. @. We focus on the 
evaluation of the normalized difference for S states, which 
we denote as 



■ AE(nS) - AE(1S) . 



(Dl) 



Important contributions to the Lamb shift as used for the data 
in Table|VIl| can be found in Tables 1—10 of Ref. O- The 
new results derived in this article for the nonlogarithmic two- 
loop term Bqo(tiS) — Bqo(1S) can now be used for an im- 
provement of the accuracy of the theoretical predictions as 
listed in TablelVTTI 

Extrapolations of the two-loop Bethe logarithms b^nS), 
and of the Aqq coefficients in Table|IJ to higher principal quan- 
tum numbers, are performed by assuming a functional form 
of the type a + b/n + c/n 2 for the correction, with a, b and 
c as constant coefficients. This functional form has recently 
been shown to be applicable to a variety of quantum electro- 
dynamic corrections for bound states, see e.g. Refs. GIGJ]. 
The same functional forms are used to extrapolate the differ- 
ence Gse(&) — Aqq, as a function of n, to higher principal 
quantum numbers [numerical results of the nonperturbative 
self-energy remainder Gse (a) can be found in Refs. iTMl30ll l. 



Although the coefficient is finite, the matrix element of ( VV) 2 
yields the logarithm 



((W) 2 )^-4^#ln[(Za)- 2 ] 
The additional contribution to Bqi is therefore 



(C3) 



SB 61 



F^(0) + (E 



(1)\ 2 



A 2 ) 



>> 



648 144 sv ; 



= -1.385 414. (C4) 



Again, the subscript vp denotes the contribution from the sub- 
set it of two-loop diagrams (Fig. |3. The new value for the 
logarithmic contribution including the vacuum polarization is 



SB 61 = 



413581 
64800 



616 

~135 
48.958 590 



N(1S) 
40 



2027 
"864" 



ln(2) - - 7T 2 ln(2) + — ln 2 (2) + £(3) 
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(C5) 



TABLE VII: Theoretical values of the normalized Lamb-shift dif- 
ference A n defined in Eq. <DH . using results obtained here [see 
Eq. <Ql. Units are kHz. 



n 


A„ 


n 


A n 


2 


187225.70(5) 


17 


281845.77(11) 


3 


235070.90(7) 


18 


282049.05(11) 


4 


254419.32(8) 


19 


282221.81(11) 


5 


264154.03(9) 


20 


282369.85(11) 


6 


269738.49(9) 


21 


282497.67(11) 


7 


273237.83(9) 


22 


282608.78(11) 


8 


275574.90(10) 


23 


282705.98(11) 


9 


277212.89(10) 


24 


282791.50(11) 


10 


278405.21(10) 


25 


282867.11(11) 


11 


279300.01(10) 


26 


282934.29(11) 


12 


279988.60(10) 


27 


282994.18(11) 


13 


280529.77(10) 


28 


283048.01(11) 


14 


280962.77(10) 


29 


283096.35(11) 


15 


281314.61(10) 


30 


283140.01(11) 


16 


281604.34(11) 


31 


283179.54(11) 



The principal theoretical uncertainty with regard to the nor- 
malized difference A„ currently originates from the unknown 
n-dependence of the two-loop coefficient BninS). An esti- 
mate for this correction may be obtained as follows. We first 



21 



map the one-loop coefficient onto an effective Dirac delta 
potential V50, with 



V f . 



a 



50 



{Zaf 



427 
384 



1 



In 2 



V 2 V. 



(D2) 



Of course, W 2 V = 47T(5 3 (r), and we may use this potential 
as an "input" for evaluation the additional one-loop correc- 
tion to the Bethe logarithm generated by the local potential. 
This leads to a correction of order a 2 (Za) 7 , with logarith- 
mic terms. The leading double-logarithmic term (correspond- 
ing to a -B72 -coefficient) is n-independent. The well-known 
n-dependence of the single logarithm, which gives rise to a 



i?7i-coefficient, may be found e.g. in Eq. (20) of Ref. H. The 
calculation leads to the estimate 



'427 1fi 

B 71 (nS) - B 71 (1S) « rr ( — - - ln(2) 



--- + — 2 + 7 + *(n)-ln(n 
4 n 4n z 



(D3) 



for the nS-lS difference of the logarithmic term. As an un- 
certainty estimate for Bn(nS) — Bn(lS), we take half the 
value of the above expression. 
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